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Abstract 

A new and conceptually simple procedure is derived for the compu- 
tation of the maximal reachability submodule of a given submodule of 
the state space of a linear discrete time system over a Noethenian ring 
R. The procedure is effective if R is effective and if kernels and inter- 
sections can be computed. The procedure is compared with a rather 
different procedure by Assan e.a. published recently. 

1 Introduction 

Let A € R^^'^, B € i?"^™ where for the moment R is just a commutative 
ring. As usual, we associate to the pair (^4, B) the linear discrete time 
control processes 

xq , xi = Axo + Buo , • • • , Xk+i = Axk + Buk , • • • (1) 

with states Xk € -R", inputs Uk € and A; € N. 

A submodule U of R^ is called ( A, -invariant if AU C U + 'miB. 
An [A, B) invariant submodule U is called reachable or reachability sub- 
module if every state in U can be reached from zero within U . The latter 
means: 

Vx € ?7 3r G N, uo,... G R"^ : 

xi = Buo, ... ,Xr = A^~^Buo -|- . . . -|- Bur-i € U and Xr = x . 



It was shown (see e.g. [|t|, Theorem 2.15]) that this rather natural definition 



is equivalent to the definition of pre-controllability submodules in |CoPe] 
which is still more commonly known but less intuitive from a control point 
of view. 

The zero-module is trivially {A, i?)-invariant and reachable. Prom the 
definitions it is clear that sums of {A, i?)-invariant or reachable submodules, 
respectively, are again (A, i?)-invariant or reachable. These facts imply that 
any submodule M of R" contains a unique maximal {A, i?)-invariant sub- 
module M* and a unique maximal reachability submodule Mg , where always 

Mg* C M*. 

Maximal reachability submodules play an important role in the solutions 
to classical control problems such as disturbance decoupling. See [ CoPe | and 



;o give only two examples. It is therefore of practical importance to 



have methods at hand for the computation of generating systems of such 



modules. In [ AsLaPel ] for the first time a finite procedure was given for 



principal ideal domains and then strongly modified in |AsLaPe2] to work for 
Noetherian rings. The latter works as follows: 

R is now supposed to be Noetherian. 

First step (precalculation): 5*0 := imB 

and for A; > 1 : Sk := imB + A{Sk-in M). 

This ascending sequence of modules stabilizes after finitely many steps 
and gives a submodule M^, which contains the image oi B. If M is 
represented as the kernel of some matrix C E i?"^^, then M^, appears 
as the 'minimal (C, A)-invariant submodule' containing the image of 
B, see e.g. |AsLaPe2| ]. 



Second step and main procedure: Wq := d M d A ^{ImB) 
and for k > 1: Wk := M^nMn A-^{Wk-i + ImB). 
Once more, this gives an ascending sequence and an interesting proof 



in [AsLaPe2| shows that its limit is actually Mq. 



Of course - and the same is valid for the new procedure to be developed in 
this note - such a procedure can be realized in a concrete computation only 
if the ring R and all the occurring operations like "n" are effective 

in the sense of []CoCuSt , p.l]. 
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2 New procedure via finite (A, S)-cyclic 
submodules 



Based on results from [BrSch| a quite different and conceptually simpler 
approach is possible. A submodule U of i?" is called (A, S)-cyclic if for 
some Uk G -R™ and Xk from (||) with Xo = one has 



U = {xk:k>0). (2) 

Thus an {A, i?)-cyclic submodule can be generated by the states of one single 
control process which begins with the zero-state. 

It is shown in [BrSch] that (^4, i?)-cyclic submodules are reachability 



submodules and that finitely generated reachability submodules are even 
finite {A, S)-cyclic. The latter means that in addition to (|2|) one has 3:^ = 
for A; > d and some d €N. 

The point is now that finite {A, i3)-cyclic submodules can be determined 
via the kernel of [yE - A, -B] in If for / G g G Rly]"" 

one has {yE — A)f = Bg, then the coefficient vectors of / generate a fi- 
nite {A, i?)-cyclic submodule and every finite {A, i3)-cyclic submodule U = 
{xi Xfi, 0, . . .) leads to a kernel element [ ^ ] with / = x\y'^~^ + . . .+ x^ 
and g = u^y'^ + . . . + u^. Note that x^+i = -^dXd + Bud = 0. More details 



can be found in | BrSch | 



For any / = xiy'^ ^ + . . . + € let Uf := {xi, ... ,Xd). 

Of course, Uj is contained in a given submodule M if and only if the co- 
efficient vectors of / are from M. Let vr be the projection of = 
© R[y]^ onto the first n components and let 

M:=KeT[yE-A,-B]n (M[y] x . (3) 

Here M[y] is the submodule of i?[y]" of those polynomial vectors which have 
all their coefficient vectors from M. 

One arrives now at the following results: 

Observation, (i) For every h £ the submodule U^(^h) is a reachability 
submodule of M (true for any R). 

(ii) Let R be Noetherian. For every reachability submodule U of M 
there is h e M such that U = C/^(/i). 



3 



Proposition. Let hi, . . . ,hs generate A4 as an R[y]-module, then the fam- 
ily of coefficient vectors of ^{hi), . . . ,vr(/is) generates Mq. 



Proof of Observation, (i): By construction UT^(^h) is finite (A, -B)-cyclic 



and thus by Proposition 1.5 in | BrScb | a reachability submodule 



(ii): Since R is Noetherian, U is finitely generated and reachable. By Propo- 



sition 1.7 in | BrSch | this implies that U is finite {A, i3)-cyclic. The foregoing 



discussion shows how to construct the desired h £ M. □ 

Proof of Proposition. Let fi = 7r(/ii), ... , fs = T^{hs) and M = Yl!i=i ^/r 
We have to show M = Mq . Mq is the sum of all reachability submodules of 
M. Since R is Noetherian, all reachability submodules U oi M are finitely 
generated. By part (ii) of the Observation such modules U can be repre- 
sented asU = U^(^h) with some h G ^A. Since h = rihi + . . .+rshs with some 
ri, . . . , € we obtain U M for an arbitrary reachability submodule 
U of M and thus Mq* C M. 

The converse inclusion comes from the fact that by part (i) of the Ob- 
servation Uf. is a reachability submodule of M and therefore contained in 
Mq* for 1 < z < s. The latter imphes: M C Mg*. □ 

One main advantage of the approach via (^) is that one can (for appro- 
priate rings R) compute the kernel of [yE—A, —B] once for all independently 
of M. This gives us as a first result a module which is of use not only for de- 
termining Mq, see e.g. |BrSch|. In order to determine Mq for some specific 



M it remains to calculate an intersection of two modules and after that one 
merely truncates the results and extracts the coefficient vectors.. Explicit 
calculation is - of course - only possible over an effective Noetherian ring 
with an effective method to determine the kernel and intersection in (^). 
Examples of such rings are Z, Q [ti, . . . , F [ii, i„] where F is a finite field. 
The determination of Ker [yE — A, —B] can then be done with the help of 
Grobner basis calculations as indicated in |BrSch|. A standard technique 



also via Grobner bases for the computation of the intersections of modules 



is (e.g.) described in [AdLou]. In both cases any generating system would 
do as well. Several current software packages for symbolic computation can 
be utilized to perform explicit calculations. 

A sound comparison of the different procedures for the computation of 
maximal reachability submodules requires a detailed investigation of their 
complexities. This remains as a future task. 
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The following two examples are over <Q[t] and Q[t,u;]. Computations 
have been done combining the well-known packages Macaulay2 and MapleV 
Release 5.1 



Examples 



(A) Let A 
of |AsLaPe2|] . 
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To determine Mq we first obtain 
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as in Example 1 



Ker[y£' - A, -B] = im 



This leads to = hR[y] with h = ^[t, —t, —t,t, —y], which in turn leads 
to with / = TT{h) = ^[t, —t, —t]. There is only one coefficient vector to be 
extracted from / (viewed as a polynomial vector in the variable y). Therefore 
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fR. By |AsLaPe2| ] we know M* 
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R and 



the final result is: 
thus Mq* C m* . 

This example is interesting also since here the classical Wonham-algorithm 
to determine M* does not converge and up to now no general finite proce- 
dure is known. For principal ideal domains, however, a procedure has been 
developed in [ AsLaLoP^ . 



(B) In the second example we start with matrices from [A.sLaPe2|, Exam- 
ple 4.3, where a system with two incommensurable delays is investigated. 

Let 





' 





1 




't 


0" 




"1 





A = 


4 

W 


t 





, B = 








, and M = im 





w 






t 


1 







1 







1 



Here Macaulay2 computes 



Ker [yE - A, -B] = im 
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-wH + f^) - t^y - ty"^ + 
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which leads to A4 = hR[y] with 



h= '[t^-ty, 



-wH, 



-w^ + ty-y^ 



-wh + t'')y] 



Now vr(/i) = Xiy + X2 where xi = *[— 1,0, 0] and *X2 = [t"^ , —w^t, —w^t] 
and according to the Proposition we obtain as final result: Mq = {xi,X2) 
(compare with i?2 ™ [A.sLaPe2, 4.3]). Note that by the new procedure we 
automatically get Mq represented as an (A, S)-cyclic subspace. In more 
complex examples one obtains Mq as a sum of {A, i?)-cyclic modules. For 
reasons of space I do not give an example for this. 
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